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To fulfill the needs of neutron capture reaction cross-section measurement in keV energy scale in the field 
of nuclear astrophysics and advanced nuclear energy system development, the 47 BaF2 Gamma-ray Total Ab- 
sorption Facility (GTAF) developed by the Key Laboratory of Nuclear Data of the China Institute of Atomic 
Energy (CIAE) has been transplanted and installed at the Back-streaming White Neutron Source (Back-n) of 
the China Spallation Neutron Source (CSNS) in 2019. A series of results has been achieved and published 
based on the GTAF since then, and due to which the needs of reducing backgrounds are becoming increasingly 
urgent. In order to understand the origins of backgrounds and to optimize the facilities, a detailed simulation 
program using GEANT4 toolkits was established and presented in this paper. To demonstrate the availability of 
the proved codes, several practical examples of assisting the process of experimental data and helping verify the 


optimization proposition are also shown in this paper. 


Keywords: Gamma-ray Total Absorption Facility, White Neutron Source, Neutron Capture Cross Section, Monte Carlo sim- 


ulation, GEANT4, Geometry optimization. 


I. INTRODUCTION 


Due to its large cover angle and high detection effi- 
ciency [1], the 4r BaF, Gamma-ray Total Absorption Facility 
(GTAP), as shown in Fig.1, is designed to meet the needs of 
neutron capture cross-section measurement under keV energy 
scale neutron beams in the topics of nuclear astrophysics and 
advanced reactor design [2—7].It has been transplanted to be 
setup at Back-n of the CSNS in 2019 [7, 8]. 


Fig. 1. GTAF Detector Array and Associated Facilities Installed in 
the Hall-2 of Back-n CSNS[2] 


In order to assist the analysis of experimental data, a set 
of detailed and reliable Monte Carlo simulation codes is pro- 
grammed using GEANT4 toolkits [9] as described in the Sec- 
tion II and Section III. A verification by standard library or 
experimental data is shown in the Section IV and Section V 
and based on which, several practical examples are presented 
in the Section VI, including the demonstrations of availabili- 
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ties in assisting the process of experimental data and in verify- 
ing the geometric optimization variants to solve the problems 
of backgrounds. 


Ii. BASIS OF FACILITY 
A. Time-of-Flight Method 


The Time of Flight (ToF) method is a commonly used 
method in measuring particles [10, 11]. It relies on the prin- 
ciple that the time it takes for a neutron to travel a known 
distance is inversely proportional to its energy, which could 
be theoretical calculated in by Equation (1) [10, 11]. 


. 72.3% L 
al By, 


where t refers to the flight time, L to the flight distance and 
En to the primary neutron energy. 

The measurement of neutron flight time is designed in 
grand accuracy [12] using specific timing hardware and soft- 
ware system [13] since it is crucial to determine the neutron’s 
energy and to reconstruct the spectrum at the GTAF [14-16]. 


(1) 


B. Miultiplicities of (n,y) Reactions 


The Multiplicity of reactions is defined as the number of 
volumes that particles have bypassed with inelastic reactions 
before being totally absorbed or escaping the sensible crystal 
array as shown in Fig.2. 

It plays a key role in benchmarking valuable information 
about reaction channels and underlying physics process [17, 
18], such as elastic scattering, inelastic scattering, radioactive 
capture, etc. since each event has a distinctive multiplicity 
signature. 
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Fig. 2. Schema of Multiple Reactions in one Event occured in Dif- 
ferent Crystals 


C. Pilled-up Energy of Event Cascades 


The de-excitation principle of isotopes in GTAF is shown 
in Fig.3. 
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(a) Principle of Isotopes De-excitation[5] 
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(b) Principle of Pilled-up Energy[2] 


Fig. 3. Principles of Event Reconstruction. 


In GTAF, to seek neutron capture reactions is of the most 
importance since it is one of the key and interest data [2, 19]. 
To distinguish the neutron capture reaction, one of the most 


54 


55 


56 


a 


7 


58 


59 


68 
69 
70 


T 


72 
73 
74 
75 


76 


77 


78 
79 
80 
8 


82 
83 
84 
85 
86 
87 
88 
89 
90 
9 


92 
93 
94 
95 
96 
97 
98 


99 


practical ways is to find the value of pilled-up released 100 
gamma-ray energy Eex since no matter how many reaction 101 
channels has been experienced, the Eex remains constant if 102 building high level accurate geometric volumes, ensuring that 


all data can be restored ideally, as shown in Equation (2). 


Ex = En + Q (2) 
where the En refers to the neutron energy and the Q to the 
reaction Q value. 


Il. MONTE CARLO SIMULATION 


A. General Idea 


As discussed in the Section I, a reliable Monte Carlo Sim- 
ulation is needed to be established in order to fulfill the needs 
of amelioration of facilities and to help of analysis the exper- 
imental data. 

The reliability of Monte Carlo simulation depends on the 
details of reconstruction in variety parts, i.e. 1) Detailed ge- 
ometry reconstruction; 2) Accurate physics configurations; 
3) Reasonable Calibration and Neutron Beam sources; 4) 
Capable event reconstruction algorithm and 5) Logical data 
restoration design. 

The GEANT4 simulation toolkits package [9] is chosen as 
it has been widely used and verified in nuclear physics and 
high energy physics with strong abilities of extensive physics 
configurations and mutual geometric reconstruction methods. 
The kernel version of GEANT4 in use in the simulation of this 
paper is 11.1.2. The general working flow of the simulation 
is shown in Fig.4. 


B. Geometry Reconstruction 


To adapt to different needs of geometry simulation under 
certain experimental conditions, different Boolean variables, 
as shown in Table 1, are offered to users as switches of geom- 
etry construction. 

The geometry is reconstructed to be with the most reason- 
able details. Apart from the mechanical fabrication errors, 
the geometric parameters and related materials of the facili- 
ties are set as same as the ones measured directly from the 
actual arrangements [20] in the Back-n of CSNS. In addition 
to using the Constructed Solid Geometry (CSG) method or 
the CSG liked methods embedded in GEANT4 toolkits, the 
subassemblies of the facilities are also reconstructed using 
CADMesh method [21] as a back-up and agile development 
option. 

The CADMesh is a valuable tool in reconstructing detec- 
tor constructions in GEANT4 simulations. It allows import- 
ing complex geometries created in Computer-Aided Design 
(CAD) software, with a support of various common ACSII 
format files, into the GEANT4 simulation program directly. 

Although both CSG and CADMesh methods are based 
on Computer Graphics geometric logics, more preset ba- 
sic graphics and logical calculation operations are provided 
by commercial CAD softwares when constructing elements 
which offer the CADMesh method possibilities of rapid 
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Fig. 4. General Data Flow of Simulation Codes 


103 the simulation of detector’s physical characteristics could ful- 
1o4 fill the crucial needs for obtaining results in particle tracking. 
105 It is especially advantageous that the CADMesh method can 
1% easily excel at handling complex 3D shapes and curved sur- 
107 faces when simulating detectors, such as for the GTAF series 
108 detectors which contain quite a lot intricate or irregular geo- 
109 metric elements. 


110 According to the topologic definition of different fields in 
11 ACSI format CAD files created by FreeCAD [22] in this 
112 paper, such as vertex positions, normals, mappings, etc., 
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as shown in Table 2 [22], an interface program is used to 
read and translate parameters to make the core program of 
GEANT4 completing the corresponding geometric construc- 
tion whereas the corresponding materials are defined subse- 
quently by the same way as CSG method in GEANT4. 


The two mentioned methods are designed to be able to 
switch between each other via a Boolean variable as shown 
in Table 1. In addition, for certain elements whose local geo- 
metric effect is no need to be detailed considered, the related 
parameters are calculated and set by an equal-volume factor. 


Table 1. Boolean Variables for Construction under Different Simulations 
Boolean Variables 


Tier N1 Tier N2 Tier N3 Default 
Calibration / / false 
Beam Test / / false 

PreConstruction_Wall / false 

PreConstruction_GTAF / false 

PreConstruction_Pipe / false 

PreConstruction Switch PreConstruction_Absorber / false 
PreConstruction_ShellSup / false 
PreConstruction_SampleSup / false 
PreConstruction_Etagere / false 

All Vacuum / false 

General Construct WallConstruct / true 
BackEnd Trap / true 

det_Shell / true 

det_Support / true 

Detector Construction det_Shelf / tue 
det_Base / true 

nPartial Crystal Construct / true 

det_Cluster / true 

/ true 

Sample Construct Gee Le 

Sample Construction Sample_Select HBE 
/ true 


SampleSup_Al false 
SampleSup_Teflon true 
SampleSup_Select false 


Sample Support Construct 


General Pipe / true 

Pipeline Construction Penetrate Pipe / true 
Middle Pipe / true 

Absorber Construction Absorber Construct Ab_Select true 
Process Model Mode / false 

Process Analysis Physics Model Mode / true 
Simple Canal I true 


Table 2. Definition of Fields of Geometry Parameters in ASCII Format Files 


Field Meaning Description 
v Coordinates of vertices Definition of a vertex through coordinate x-y-z data in each line. 
vt Coordinates of vertex texture Definition of a vertex texture through coordinate x-y data. 


Definition of normal (number of normals is determined by 
the intersection of each vertex and face) 
In Computer Graphics, mesh is used as the definition of faces. 
f Face Every three points on different lines at least with three index values: 
vertex, vertex texture and normal could define a face. 


vn List of vertex normal 


o Objects 
Groups 
s Smoothing group 
123 For instance, the bellows type BP300 can be considered as a 129 C. Physics Models 


124 tube with a volume equivalent coefficient of 1.2 while quick 
125 release flanges type KF100 would be of 1.47. The geometric 


12 simulation and some of the typical subassemblies are shown 1% In order to simulate the detector, it is crucial to construct a 
28 in Fig.5. 131 reasonable physical model. Considering the simulation under 


132 NO matter energy calibration mode or neutron beam mode, 
133 More interests are concentrated on the response in low en- 
134 ergy range, the QGSP_BIC_HP preset physics package [23] 
135 is used as a basic physics model package, which contains a 
e series of physics reference including the low-energy reaction, 
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(a) Mass Plan of the Hall-2 Geometry 
using CSG+CAD Method 
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Fig. 5. Typical Construction of Facilities 


decay, elastic scattering and inelastic processes that can meet 
the preliminary needs of simulation with results shown in 
Section V in this paper. The "HP" type physical process pack- 
age is used here since the character "HP" refers to High Pre- 
cision physics models [23] in the context of GEANT4 which 
provide more accurate and detailed simulations of particle in- 
teractions covered with a broader energy range for particles 
and allowed for more comprehensive simulations of various 
physics processes. The preliminary applied physics processes 
and models are shown in Table 3 and there will be continued 
to be refined in subsequent work due to different needs of 
simulation. 


Table 3. Applied Physics Processes and Models 


; Mode 
Physics Model Calib Neutron 
QGSP_ BIC_HP a a 
EMV_option4 a E 
DecayPhysics E a 
BiasedRDPhysics E 
HardronElasticPhysicsHP E E 
IonElasticPhysics a 
IonPhysics E 
GammaNuclearPhysics E 
GammaNuclearPhysicsLEND m 
(option) 
NeutronHPPhysics E 
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D. Primary Sources 


Two types of sources, i.e., 1) Calibration sources recom- 
mended by the standard library Evaluated Nuclear Data File 
(ENDF) [24, 25]: ®°Co, 137Cs, and ??Na; 2) Neutron beams: 
including a) Neutron beams output from Back-n and b) 4.9eV 
mono-energetic neutron beams, are reconstructed using spe- 
cific macro files, in which contains a matrix of spectral and 
spatial parameters and related normalized weighting coeffi- 
cients to provide the required information for the simulation 
of primary sources. 


1. Calibration Sources 


Three Calibration Sources recommended by ENDF library, 
60Co, 187Cs and 22Na, are reconstructed in the simulation 
with the same geometric parameters as used in real experi- 
ments [26]. The dimensions are set as $32 x 4mm, ¢32 x 4mm 
and 25 x 25mm respectively. 


2. Neutron Sources 


Two neutron beam sources, mono-energetic neutron beams 
and spectral neutron beams are simulated. Due to the time- 
resolution limit of hardware in real experiments, Time-of- 
Flight (ToF) spectrum of beams with initial energy over 1 
MeV cannot be well resolved. Therefore, on the first stage 
of simulation, an upper energy limit of 1 MeV has been set to 
both of the options [26]. The matrix of neutron energy spec- 
trum and initial momentum spectrum is written in a macro 
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1. The mono-energetic neutron beams are simulated with 
parameters of a 4.9 eV mono-energy and a spatial 
Beam spot obtained from the CMOS experiment data 
at the Back-n of CSNS. It is used for obtaining a clear 
image of the largest resonant cross-section of the stan- 
dard 197Au sample which leads to verify the reliability 
of the codes and to be used to help calculate the theo- 
retical efficiency or other required information. 


. The spectral neutron beams are simulated with the 
same spectral [27, 28] and spatial[29, 30] characteris- 
tics of which from the Back-n. Itis used for the anal- 
ysis of backgrounds and for the calculation of the the- 
oretical neutron capture cross-section of samples. The 
simulated beam spot is shown in Fig.6. 


3. Pseudo-random Number Generator 


The statistical properties of Pseudo-Random Number Gen- 
erator (PRNG) have a great impact on the reliability of Monte 
Carlo simulation results. There are several popular candidate 
PRNG s for nuclear physics, as the James Random, Mersenne 
Twister and Ranlux64. The MTWistEngine pseudo-random 
number generator[31] is chosen as the PRNG used in the sim- 
ulation due to two main reasons, as: 


1. a big enough pool of valid pseudo-random number of 
219937 — 1 can be generated at one single operation 
which might support the needs of an upper limit of 
about 2 x 10° events in each run in GEANT4 toolk- 
its and might also fulfill the potential needs for further 
study using the accumulated simulation data; 


2. ahigh reliability since it has been passed almost all the 
rigorous random number tests referred to the analysis 


thesis in reference [32-34]. 


In order to get the system beeper via an I/O data channel 
as initial seeds feeding the pseudo-random generator, a time 
seed interface program is designed, programmed and linked 
to the EventAction class of the simulation program. 


E. Simulation Run and Action Classes 


Under no matter the Calibration mode or the Neutron Beam 
mode, the simulation can be started by emitting the primary 
particles by calling the detail parameters of the initial parti- 
cles in the macro file (including the initial energy, initial mo- 
mentum, and initial position) and the specific operation mode 
defined in the RunMessenger interface, via the General Parti- 
cle Sources (GPS) function embedded in GEANT4. 

After emission, the behavior of particles depends on the 
different physical models and definitions of the physical pro- 
cesses in Section III C while travelling through various spatial 
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volumes in the geometry described in Section IIIB and per- 
forms to cutoff or truncate the information according to the 
demands in the various action classes. 

According to the design logic of action classes built-in 
GEANT4, each Run consists of the number of Events de- 
fined in the macro file, and in each Event, according to the 
setting of the particle step length and physical process, the 
corresponding Step is included, and the required information 
can be filtered through the action classes of Stack and Track. 

At the end of each Run, i.e. after the last Event, the cor- 
responding data and information will be saved and output ac- 
cording to the required format preset in RunMessenger. 


IV. PRIMARY DATA ANALYSIS PROGRAM 


In order to meet the basic pre-processing requirements of 
the simulation data, a general data processing program with a 
set of Qt-based GUI interface is designed and tested. 


A. GUI Interface 


In order to facilitate the implementation of commonly 
used data preliminary processing functions, a visual human- 
machine GUI interface through Qt version 5.9.7 is designed, 
as shown in Fig.7. 

With PythonQt and PyRoot, the commonly used functions 
are realized by transporting the data flow via various inter- 
face, including the import and export of files in various basic 
formats, switching between neutron beam mode and calibra- 
tion mode, display and fitting functions of energy spectrum 
and ToF spectrum. 

The visualization program uses the control of Qt and re- 
alizes the data interface required to call the above functions, 
and realizes the basic data preliminary processing. 


B. Reconstruction Algorithm of Event Cascades 


The simulation can be done by each event, or each run as 
required. Consistent with the process of experimental data 
processing, two general and basic event reconstruction algo- 
rithm subprograms in the data processing program are de- 
signed: energy reconstruction and position reconstruction. 

The reconstruction is relatively simpler in the simulation 
since it is possible to retrieve the target data directly from 
GEANT4 built-in functions. Particles are transported and 
tracked via the functions of Action Classes in the framework 
of GEANT4 toolkits until they are absorbed in certain vol- 
umes or escape the set cut-off areas. 

The essential value, as the deposition energy, the time-of- 
flight, the multiplicities, the geometric volumes, the material, 
the reaction channel, and other relevant information can be 
taken from and recorded after each step or event. Meanwhile, 
the corresponding data of each reaction channel can be dis- 
tinguished by calling the physics model or physics process of 
each track. 
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Fig. 6. Simulation of Neutron Beam Spot at Back-n of CSNS 
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Fig. 7. GUI Interface of GTAF Simulation Pre-processing Program 


1. Energy Reconstruction 


The deposited energy could be traced in each step of 
particles transportation with help of build-in algorithms of 
GEANT4 which makes one of the privileges of processing 
simulations data compared to the counterpart of experiment 
data is that the pilled-up energy peaks and event cascades re- 
constructions are much easier to be done. 

The particles are designed to be transported via the func- 
tions of action class in each Step until they are absorbed in 
certain volumes or escape the preset cut-off areas. While the 
simulation runs, each particle will be labeled by each gener- 
ation, the associated information of particles would be trans- 
ferred via a user-set personalized function to the TrackingAc- 
tion class in order to get further processed. 

Considering that the simulation results of interest are those 
with responses in the detector array, the geometric volume 
information is of more importance. Therefore, a label each 
volume number (the CopyID) are recorded at the same time. 

For situations where further processing of energy deposi- 
tion is not required, the above-extracted data can be directly 
transferred from the SteppingAction class to the EventAction 
class so as to be accumulated and stored directly. And after 
which, the data would be finally passed to the Analyzer for 
saving in a certain format according to the designed dataflow 
shown in Fig.4. 

On the other hand, where there is a need for further pro- 
cessing, a TrackingAction class retains interface functions 
for users to filter the specific required deposition energy is 
designed. According to the needs of the user’s simulation 
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objectives, data can be transferred to the Analyzer after pre- 
liminary processed in the TrackingAction class, where more 
preprogrammed reorganization tools, such as different reac- 
tion channels, different multiplicities, and different areas to 
store and reconstruct the required energy spectrum, are set as 
to be described in the Section IV C. 

The original simulation data output from the dataflow 
above are saved by separate detector crystals. The recon- 
structed energy spectrum can be output directly in the form of 
divided crystals, or according to the user’s needs, a deposited 
energy spectrum reconstruction output in the unit of the total 
detector can be achieved in the form of the sum of the total 
deposition data of each crystal in the same Event. 

The above functions can be implemented in the RunMes- 
senger by adjusting relative Boolean variables, or through 
GUI tabs. Thus, after each Run, the energy spectrum of par- 
ticles could be reconstructed. 


2. Time-of-Flight Spectrum Simulation 


In order to verify the results with experimental data, the 
flight time of particles is designed to be recorded in the simu- 
lation. 

The very start time point TO is preset in the EventAction 
class at the beginning of Event when the primary particles 
begin to emit in every loop. When triggered in each Step 
(under SteppingAction mode) or in each Sensitive Detector 
zone (under SensitiveDetector mode), the corresponding time 
is recorded and saved in a tuple or histogram predeclared in 
the RunAction class. Note that the recorded time mentioned 
above is a Global Time in the entire Event since TO is the 
beginning of each Event. Right at the end of each Event, 
corresponding time data are recorded in different tuples or 
trees in the ROOT files through the pre-selection conditions in 
action classes of each Step, Track, Stack or Event. Therefore, 
a ToF Spectrum can be generated at the end of the whole Run, 
e.g., the end of the final Event. 

In addition, similar to processing of experimental data, the 
corresponding simulated energy spectrum (E-ToF) can be cal- 
culated using the flight length of particles and the simulated 
ToF spectrum through the mentioned Equation 1. The flight 
length of particles is obtained by adding the flight distance in 
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each Step. It is calculated in geometric simulation program 
and transferred to the analysis functions by two following 
ways which can be switched into each other by users through 
a Boolean variable. 


1. Calculated by extracting the geometric length of the 
corresponding passing geometric elements in the De- 
tector Construction source file, and the very geometric 
parameters would be transferred to the Analyzer source 
file through the transfer function; 


2. Calculated directly in Step action class through the 
built-in variable function of GEANT4 toolkit, where- 
after the step length would be passed to the function in 
the Event action class in order to store and generate the 
E-ToF spectrum directly. 


3. Position Reconstruction 


Similar to the reconstruction of ToF spectrum, the position 
information (3-Dimensions vector tuple) of each step can be 
traced and recorded while the deposited energy (the differ- 
ence between the pre-step and post-step energy) in dedicated 
geometric volumes is not equaled to zero in the SteppingAc- 
tion class. 


C. Reorganization Tools 
1. Multiplicities 


During the experiment installation setup stage, as discussed 
in the Section II B, a proper design should be considered for 
processing and identification the multiplicities. In GTAF, 
an electronic circuitry based on NOT gate circuit of nuclear 
electronics technology is realized with several key parame- 
ters well preset including the energy and timing thresholds to 
identify particle interactions [14]. 

In the simulation program, a similar but of more precision 
and practical method has also been applied. Technically, the 
multiplicities of each Event are counted by the number of dif- 
ferent CopyIDs of geometric volumes that the deposit energy 
is not equal to zero before the particle is fully absorbed or es- 
capes from the sensible arrays since there is a unique CopyID 
tagged to each of the geometric volumes reconstructed in the 
simulation codes. 


2. Reaction Channels 


In experimental data processing, distinguishing data from 
different reaction channels is the core algorithm of data pro- 
cessing, which is realized through different gates in order to 
help understand the experimental data and phenomena. 

In simulation, the physical process occurred at each step 
could be traced before or after each step. In order to avoid 
the null Pointer error in C++ coding environment, apart from 


3 


foo} 


5 
386 
387 
388 
389 
390 
391 
392 
393 
394 
395 
396 
397 
398 


399 


400 


401 


402 
403 
404 
405 
406 
407 
408 


4 


S 


9 


410 
411 


412 


413 


414 


415 


416 
417 
418 


419 


420 


421 


422 
423 
424 
425 
426 
427 
428 
429 
430 
43 


432 


the usual protection by a judgement function, the post step 
physics model filter is used in this program. A string value, 
with a Pronouns preset in GEANT4 or by users-set to the ded- 
icated physics model or physics process, due to the Boolean 
values switch that has been chosen in the RunManager, would 
be returned. After transferred to the Stack and Track classes 
of the simulation program, the string values of the relevant 
physics processes or physics models are passed to the Ana- 
lyzer and stored in the corresponding tuple or other format in 
a ROOT file. 

After the simulation is completed, the value of the physics 
processes or physics models could be called in Primary Data 
Analysis Program, and the ToF spectrum or energy spectrum 
involved in different reaction channels can be classified and 
plotted. 


D. Spectrum Broadening and Semi-automatic Peak Finding 
1. Broadening of the Energy Spectrum 


Since GEANT4 cannot simulate the nuclear electronic ef- 
fect in the preset physical process, the ratio of electronic re- 
sponse is obviously 100%. Thus, the simulated data need to 
be broadened before supporting the experimental data analy- 
sis. 

In the preliminary analysis program, Gaussian’s function is 
used as the broadening algorithm. The specific process of the 
broadening algorithm is as follows: 


1. Determine the total normalized bin number and the cor- 
responding coordinate value of the corresponding spec- 
trum (or the corresponding segment of the spectrum); 


2. Determine an energy resolution that is set according to 
the experiment or set by the user; 


3. Determine the width of the error limit; 


. Determine the constant of Gaussian broadening: by en- 
suring that the integral of Gaussian broadening with 
the above parameters is the same as the original count 
value; 


5. Plotting and recording the parameters. 


2. Semi-automatic Fitting and Peak Finding 


Spectrum fitting and peak finding are generally performed 
in the range of interest when processing data. The semi- 
automatic spectrum fitting, and peak finding process can be 
implemented in the primary data analysis program. 

The algorithm of spectrum fitting, and peak finding is sim- 
ilar to the counterpart of energy spectrum broadening, which 
is also achieved by fitting with Gaussian function with a series 
of basic initial parameters, including the approximate region 
of the peak position, basic fitting adjustment parameters, etc. 
that could be modified and fitted directly in the GUI interface 
mentioned in Section IV A. 
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The peak position and the final coefficient of the fitting it- 
eration will be displayed in the display area of the GUI inter- 
face or be printed in the Terminal, which will be stored and 
be used for subsequent data analysis. 


V. VALIDATION OF RELIABILITY 
A. Responses to Calibration Source 


Three simulated calibration sources mentioned in the Sec- 
tion IID are designed to validate the reliability of the geo- 
metric simulations and the algorithm of reconstruction. The 
results are shown in Fig.8, in which the peaks of piled-up de- 
posit energy are all in good agreement with data from ENDF 
library that well demonstrates the reliability of the geometry 
and physics configurations. 


Simulation Calibration Energie Accumulation Spectrum 
TE 


Counts 


2.5 
Deposit Enérgie/MeV 


Fig. 8. Calibration Results of Typical Simulation 


In addition, preliminary processing of multiplicities and 
reunited-BaF_2 crystals event reconstruction is performed 
for the simulation data. Taking the simulation data of °°Co 
source calibration as an example, two gamma-rays with en- 
ergy of 1.17 MeV and 1.33 MeV respectively, emit sponta- 
neously. The pilled-up energy of 2.5 MeV could be a bench- 
mark to evaluate the efficiency of the detector array as dis- 
cussed in the Section III E, as shown in Fig.9. 

It could preliminarily be proved that the geometric recon- 
struction of the simulation program is effective, and the basic 
reconstruction algorithm is available. At present, the experi- 
mental data processing of GTAF is still ongoing and the con- 
trol results from the experiment sides would be published in 
consequence. 


1. Response to Neutron Capture Reactions 


Considering that there is a very large resonant neutron cap- 
ture cross-section at 4.9 eV of the isotope 197Au, which is 
several orders of magnitude larger than other cross-sections 
such as elastic scattering, the very specific monoenergetic 
neutron beam is commonly used for verifying the physics 
configurations of the simulation. 

To this end, a 4.9 eV monoenergetic neutron beam with 
the same geometric spatial distribution as the Back-n neutron 
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(b) Energy Spectrum under Different Multiplicities Filters 
Conditions 


Fig. 9. Demonstration of Multiplicities Identification for Simulated 
Experiment of ĉ°Co Sample 


source starts from the vacuum tube 72.7 m up away the sam- 
ple tray and a standard thin cylindrical 197Au Sample with a 
same geometric dimension as experiment, i.e. a thickness of 
0.2 mm and a diameter of 40 mm, is simulated. 

In this simulation, a lower energy threshold of 10‘~?) eV is 
preset for each simulated crystal unit to facilitate preliminary 
data processing. The response information of particles on the 
BaF2 crystal was recorded as described in the Section IV B 
and after which, the energy spectrum and ToF spectrum were 
output through event reconstruction. 

As shown in Fig. 10 a peak of deposited energy around 6.51 
MeV and a typical time peak of 2.478 x 10° ns are clearly 
shown in the Energy Spectrum and in the ToF Spectrum re- 
spectively, which are consistent with the standard values and 
demonstrated well the validation of the simulation codes. 


VI. PRACTICAL EXAMPLES 


A. Assistant Processing and Understanding Experimental 
Data 


1. Impact of Different Neutron Beam Energy 


In order to speed up the simulation and considering the 
characteristics of electronic devices in real experiment, the 
Back-n neutron beam energy segment below 1 MeV is often 
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Fig. 10. Simulation Results of °” Au Sample Response to 4.9 eV 
Monoenergetic Neutron Beam 


used as the input neutron beam source. However, neutron 
beams in different energy bands may have different effects on 
the background [35]. To confirm the influence of the high- 
energy band on the part of the effect of interest, four different 
initial input neutron beam sources are simulated respectively, 
with Effect-Background Ratio results shown in Table 4. 


As can be seen from the above table, although the high- 
energy segment has a certain influence on the spectral struc- 
ture, it has little effect on the key data, as the background- 
effect ratio. Therefore, when using simulation calculations 
for rough analysis, a simplified neutron source term, with a 
filter of under 1MeV energy spectrum, can be used to improve 
the computational efficiency. 


2. Discrimination of Different Reaction Channels 


The Discrimination of different reaction channels is of the 
most privileges to use the simulation codes since it could pro- 
vide an ideal panorama of all reactions occurred. A demon- 
stration of this function is shown in Fig. 11. 


By realizing the presented functions, theoretical neutron 
capture reaction detecting efficiency can be calculated. Be- 
sides, it is an important analysis tool to provide possibilities 
to help better understand the phenomena of experimental data 
and optimize the structure by reducing the background im- 
pact. 
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(b) Energy Spectrum (Different Reaction Channels) 


Fig. 11. Demonstrations of Reaction Channels Discrimination 


B. Assistant in Evaluation of Preliminary Geometric 
Optimization 


1. Theoretical Analysis of Background 


In order to support the coming upgrade of the facilities, 
the theoretical backgrounds could be analyzed with the help 
of the simulation codes while comparing with the experiment 
results obtained [36]. 

With all of the experiment facilities installed in the Hall 2 
being simulated, since part of background could be produced 
by the interactions of scattering neutrons with surroundings 
[37], according to Fig.11, a series of abnormal resonant peaks 
displayed in the ToF spectrum ranged from 8x10° ns to 
1.1x10° ns. The preceding geometric volumes and related 
materials of the abnormal data are traced by the simulation 
codes, as shown in Fig.12. 

Considering that most of the precursor origins of back- 
grounds are the concrete-made volumes (walls, ceilings or 
floors), a preliminary proposition of geometric optimization 
thus could be made as one possible proposition to reduce the 
impact of abnormal backgrounds is to isolate the sub-particles 
caused by scattering by the wall, the ceiling, or the floor, es- 
pecially those which would impact the response in the central 
area where the crystal array lied. 


2. Evaluation of Geometric Optimization Proposition 


One possible structural optimization option is to add a vac- 
uum tube in the central area (i.e., the area through the center 
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Table 4. Theoretical Effect-Background Ratio under Different Simulation Conditions 


Neutron Beam Condition 


Number of Simulation 
Energy Spectrum 


Spatial Structure 


Effect-Background Ratio 


Back-n Energy Spectrum 


N_04 (filter under 1MeV) Back-n Spatial Structure 7.26% 

N_29 Back-n Energy Spectrum Back-n Spatial Structure 7.11% 
ces ss» Which allows us to use it in assisting the analysis of experi- 

Erma . . . “a4: 
seo mental data and helping optimize the facilities. 
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(b) Source of Precursor Material Analysis 


Fig. 12. Demonstration of Preliminary Analysis of Backgrounds 


of detector array) with a ball-shaped neutron absorber out- 
side the sample tray/support, as shown in Fig.13. The simu- 
lation results for the different preliminarily geometric propo- 
sitions are summarized in Table 5, in which show that the 
addition of the central vacuum pipe and the absorber outside 
the sample tray can significantly help in reducing the influ- 
ence of the anomalous background. Obviously, the final geo- 
metric optimization plan would be decided after considering 
more details, including the effects of in-beam y rays at Back- 
n sources[38]. Several validation experiments are prepared to 
be done, and after which all the simulated and experimental 
data would be verified and analyzed in the very near future. 


VII. SUMMARY 


A Monte Carlo simulation program for the GTAF based 
on GEANT4 toolkits is established and verified in this paper 
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(a) Mass Plan with the Optimized Geometry 


(b) Detail simulation of the Optimized Geometry in Central 
Zone 


Fig. 13. Mass Plan with the Optimized Geometry 


The geometry of the entire facilities is reconstructed in 
great detail according to the as-build drawings and the actual 
layout conditions on site. Together with reasonable physics 
configurations and event reconstruction algorithms, the codes 
have been tested and validated by comparing simulated with 
experimental data of three types of calibration sources and 
two types of neutron beam sources. All of the comparison 
results show positive agreements which demonstrate the reli- 
ability of the codes created. 

Two types of typical application examples are presented 
at the end of this paper to show some commonly scenarios 
where the above codes can be applied. 

More work will be done to enforce the performance of the 
codes and more applicable scenarios will be developed to 
help data analysis and other requests by the proved simula- 
tion codes. 
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Table 5. Theoretical Effect-Background Ratio of Different Simulations 


Central Pipe 


Absorber outside 


Number of Simulation Sample Tray Effect-Background Ratio 
Material Dimension Material Dimension 
N_04 N/A N/A N/A N/A 7.26% 
N_08 Stainless Steel 304 51 51 15.87% 
N_09 ; 52 Polyethylene (30% boron carbide) 52 15.31% 
N_10 Aluminum Alloy 6061 55 55 16% 
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